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INFERENCE OF TEMPERATURE AND WATER -VAPOR STRUCTURE 
IN THE STRATOSPHERE FROM LIMB RADIANCE PROFILES 

By Thomas B. McKee 
Langley Research Center 

SUMMARY 

Techniques to infer stratospheric temperature and water-vapor mixing ratio as a 
function of pressure from limb radiance profiles have been presented. Requirements for 
accurate pointing information have been eliminated, and only the easily obtained vertical 
scan rate is needed. Example calculations of inferred temperatures and water -vapor 
mixing ratios from unperturbed and perturbed radiance profiles are given. Radiance 
profiles for the 615 cm“^ to 715 cm"* spectral region due to thermal emission of carbon 
dioxide were used for temperature inference, and radiance profiles for the 205 cm"l to 
295 cm"l spectral region due to thermal emission of water vapor were used for mixing- 
ratio inference. Pressure ranges of the inferred values are 0.3 to 100 millibars (30 to 
10 000 N/m^) for temperature and 3 to 200 millibars for wafer-vapor mixing ratio. 
Errors in inferred quantities due to radiance errors have been examined and are not 
excessive. 


INTRODUCTION 

Techniques to infer atmospheric temperature and water -vapor structure in the 
stratosphere from limb (horizon) radiance profiles have been developed and demonstrated. 
References 1 and 2 discuss inference of temperature, illustrate the different approaches 
taken to obtain temperature structure from limb radiance profiles, and present examples 
based on measured data. References 2 and 3 discuss techniques to infer water -vapor 
mixing ratio and present examples based on measured radiance profiles. 

All of the techniques previously reported require radiance as a function of tangent 
height in order to infer temperature or water-vapor mixing ratio as a function of geo- 
metric altitude. The need to measure tangent height accurately places a severe restraint 
on any system of measuring limb radiance profiles. A pointing accuracy of 1 km at the 
limb from an altitude of 500 km corresponds to an angular accuracy of 0.023 degree, 
which is difficult to achieve. Inferring atmospheric characteristics from limb radiance 
profiles would clearly be a much simpler task if the need for accurate pointing were 
eliminated. 



The purpose of this report is to derive an inference technique which does not 
require accurate pointing data. A technique of inferring temperature and water -vapor 
mixing ratio as functions of pressure from limb radiance profiles without knowledge of 
tangent height is presented. 

SYMBOLS 

Cj constant, 1.1909 x 10"^ erg-centimeter^/second-steradian 

C 2 constant, 1.4389 centimeters -degrees Kelvin 

g acceleration of gravity, kilometers/second^ 

H tangent height, kilometers 

h altitude, kilometers 

J v source function, watts/m eter^-steradian-centimeter"^ 

N radiance, watts/meter 2 -steradian 

N c calculated radiance, watts / meter^ - steradian 

N m measured radiance, watts/meter^ -steradian 

Ny spectral radiance, watts/meter2-steradian-centimeter"l 

P pressure at tangent point, millibars 

p pressure, millibars (l millibar = 102 newtons/meter2) 

R gas constant for air, 2.87 x 10^ erg/gram-degree Kelvin 

r earth radius, kilometers 

s distance along line of sight, kilometers 

T temperature, degrees Kelvin 
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T av average temperature for layer, degrees Kelvin 

t time, seconds 

y lapse rate, degrees Kelvin/kilometer 

6 nadir angle, degrees 

v wave number, centimeter 

i >i lower wave -number limit of integration, centimeter 

v u upper wave -number limit of integration, centimeter "I 

a standard deviation 

r transmittance 

Subscripts: 

a value at center of wave band 

i integer 

o value at upper boundary of iteration technique 

1,2, 3, 4 points along line of sight 

DERIVATION OF TECHNIQUE 

The technique employed to infer atmospheric temperature and constituent mixing- 
ratio structure is an iterative process applied to the equation of radiative transfer used 
to calculate limb radiance profiles. A discussion of the calculation of radiance profiles 
and the derivation of the technique for inference are presented. 

Calculation of Radiance Profiles 

An instrument (radiometer) with a small field of view is positioned outside the 
atmosphere to measure limb radiance profiles as it scans across the limb. One line of 
sight through the atmosphere is shown in figure 1. The radiance exiting the atmosphere 



along a line of sight is given by the equation of radiative transfer (see ref. 4) as 

N = -f f Jj, — ds + C Jj/ri &v (1) 

J l> 

An evaluation of equation (1) along different lines of sight yields a radiance profile. In 
reference 1 each line of sight is defined by its tangent height H. However, the point 
where the line of sight forms a right angle with the radius is also uniquely defined by the 
total atmospheric pressure P at that point. For the present discussion, a limb radi- 
ance profile is defined as the radiance from different lines of sight as a function of pres- 
sure at the tangent point. 

Two restrictions applied to equation (1) are that no boundaries (earth or clouds) 
are encountered along the line of sight and that only thermal radiation from atmospheric 
gases contributes to the source function J v . Then equation (1) becomes 

N<P) = -IJs N ‘'S dst “' (2) 


where 


Ny = 


C ^ 3 



The meteorological variables - temperature, pressure, and mixing ratio - are needed 
to evaluate equation (2). The spectral radiance Nj, is a function of temperature only. 
The transmittance is weakly dependent on temperature but strongly dependent on pres- 
sure and the amount of absorbing gas; thus, transmittance is primarily influenced by 
mixing ratio and pressure. 

The integration of equation (2) over the variables can be transformed to an integra- 
tion of that equation over pressure p. From figure 1 the distance from the tangent point 
outward is given by 


s 2 = (r + h) 2 - (r + H) 2 = 2r(h - H) + h 2 - H 2 (3) 

Differentiation of equation (3) with r and H considered constants for one line of sight 
yields 


s ds = r dli + h dli 


(4a) 
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An approximation is possible here since r » h. The range of h involved is where 
thermal emission by carbon dioxide and water vapor is significant. Because very little 
radiation originates above 70 km, h is only 1 percent or less of r; thus, the second 
term of equation 4(a) can be omitted to yield 

s ds = r dh (4b) 

If the atmosphere is considered to be an ideal gas in hydrostatic equilibrium, the hydro- 
static equation is 


dp = - 


pg 

RT 


dh 


Combining equation (4b) and equation (5) yields 


(5) 


s ds = - — — - dp (6) 

Pg 

which is integrated to 



Equations (6) and (7) are used with equation (2) to compute limb radiance profiles without 
using altitudes. A radiance profile computed for the 615-cm - l to 715-cm - ^ spectral 
region due to thermal emission of carbon dioxide is shown in figure 2. Computations 
were made by using the 1962 U.S. Standard Atmosphere (ref. 5) and a carbon dioxide 
mixing ratio of 314 parts per million. Equation (2) was evaluated, as described in ref- 
erence 6, with transmittance data from reference 7. Doppler broadening was accounted 
for. In figure 2 the values computed by using pressure are compared with those com- 
puted by using altitude (ref. 1). 

A radiance profile computed by using pressure for the 205-cm - ^ to 295-cm - ^ 
spectral band due to thermal emission of water vapor is shown in figure 3. Computations 
were made by using the 1962 U.S. Standard Atmosphere with the mixing ratios shown in 
figure 3. Transmittance data were obtained from reference 8. 
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Inference of Temperature 


Inference of temperature requires three assumptions: 

(1) Transmittance model 

(2) Mixing ratio as a function of pressure 

(3) Absence of clouds 

Examples given in the present discussion are for the 615-cm"-*- to 715-cm - ^ spectral 
region of the carbon dioxide emission, but any suitable gas and spectral region could be 
used. 

The iterative technique to obtain temperature begins at small radiance values and 
proceeds to larger values. At the first level to be inferred, the radiance is calculated as 

r v u r T i 

N c P 0 ) = -\ N„(T)drdi, (8) 

Jv l Jl 


where the atmosphere above P 0 is assumed to be hydrostatic and have a constant lapse 
rate with height such that T at any p is given by 

-Ry 

T = T 0 ^) g (9) 

Initial values of T 0 , P Q , and y are assumed based on correlation of these parameters 
with radiances calculated from the climatological atmospheres in reference 9. From a 
measured radiance profile, the radiance error between calculated and measured values 
is 


AN(P 0 ) = N m (P Q ) -N C (P 0 ) (10) 

At this point a AT, which is a function of AN(p 0 ), is calculated, as in reference 1, by 

§» m (P 0 ) - Nc(p 0 )]t 0 2 


AT = 


/ \ C Vu f Tl 

N^,a(T 0 C2^a\ l * & 

u Vi 


( 11 ) 


where the limits r = 1 to T]_ account for the transmittance through the atmosphere. 
When T 0 is changed by this AT, the new values of T above the P Q level are 
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computed from equation (9). The values of P 0 and y remain unchanged. A new value 
of N c (P 0 ) is computed, and this procedure is repeated until AN(p 0 ) is made small. 

The calculation proceeds to the next larger radiance value. To calculate a radi- 
ance value for comparison with the measured radiance, the location of the line of sight 
in the atmosphere must be known. An examination of the radiance measurement, the 
geometry, and the hydrostatic condition enables specification of the pressure at the tan- 
gent point of the line of sight without accurate pointing information. A limb radiance 
profile is generated as an instrument scans across the limb as a function of time. The 
geometry yields a relationship between tangent height and nadir angle as 

r + H = (r + h) sin 9 

For a limb scan, the radius is constant and observer altitude is constant so that the 
time derivative is 


^ = (r + h) cos 9 ^ (12) 

The value of r + h is determined from tracking data, and cos 9 can be predicted to an 
accuracy of 2 percent or better by knowledge of the spectral region of the measurement 
for observer altitudes of 500 km or higher. Then the vertical scan rate d0/dt, which is 
an easily obtained parameter, must be measured for the experiment in order that dH/dt 
can be determined. Now, the hydrostatic equation can be expressed in terms of the vari- 
ables which define the tangent point as 


dP = - 


Pg 

RT 


dH 


For small increments of dP and dH, the upper boundary is Pj_j, the lower is P i? 
P* + P 4 1 

and the average is — — • 1 ~- . Then from the hydrostatic equation 

a 


p _ p ^ RT av + g^H 

i - i-1 2 RT av - g AH 

is obtained. The change in tangent height is AH and is expressed by (dH/dt)At which 
leads to 

2RT av + g(dH/dt)At 

1 = i_1 2RT av " g(dH/dt)At ( 13 ) 
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Equation (13) provides the specification of the pressure of the tangent point needed to 
calculate a radiance and insures that the layer is in hydrostatic balance. The AH used 
must be kept small so that equation (13), still a differential form of the hydrostatic equa- 
tion, is valid. 

Equation (13) is used with Pj_j = P 0 and T av = T 0 to estimate the new pres- 
sure P^. A radiance is next computed for the ith layer and the radiance error is 
expressed as 

ANi(Pi) = N mji (P i ) - N Cji (Pi) (14) 

From ANi, a ATi is calculated, as in reference 1, to be 

kW - N e (Pi)]Ti 2 

ATi= — ■ (15 > 

K i-,a( T l) c 2"a \ \ dT d >' 

Jt 3, 

where the limits and T3 are used; these limits account only for the transmittance 
of the new layer. Equations (14), (15), and (13) are used in an iterative operation until 
ANj is small and the new layer is in hydrostatic balance. The calculation then proceeds 
to the next radiance value and, thus, works down through the atmosphere layer by layer. 

Inference of Water -Vapor Miking Ratio 
Inference of water-vapor mixing ratio requires four assumptions: 

(1) Transmittance model 

(2) Temperature as a function of pressure 

(3) Central pressure of line of sight 

(4) Absence of clouds 

Assumptions (2) and (3) are met by optically alining the radiometer used for measuring 
water-vapor emission with the radiometer used to gather data for temperature inference. 
Water-vapor mixing ratio is assumed constant above the highest level for which radiance 
data are available. A radiance is calculated from equation (2) and a radiance error 
formed for the first measured radiance as 

* N ( p o)=N m (P 0 )-N c (P 0 ) 

The mixing ratio is changed and AN (P 0 ) determined until AN(P 0 ) is made small. 
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The next radiance point is selected and a radiance calculated for the new line of 
sight. Equation (2) is expanded to give 

N( p i) = - J J 3 Nj, dr dv - J J 2 dr dv - J J * N v dr dv (16) 

The integral from T 3 to T 2 gives a contribution from the new layer which was not . 
present in the line of sight designated by P 0 . Mixing ratio in the segment T 3 to t% 
is the only unknown needed to evaluate equation (16). A mixing ratio equal to the one in 
the layer immediately above is initially chosen to compute an N c (PjJ which is used to 
form the radiance error of 


AN(P i )=N m (p i )-N c (P i ) 

Mixing ratio is changed until AN(Pj) is made small, and the inference moves downward 
layer by layer. The method to infer water-vapor mixing ratio is the same in principle 
as that reported in reference 3. 


DISCUSSION 

I 

The technique derived herein is ideally suited to a satellite experiment for mapping 
temperature and water vapor in the stratosphere over the entire earth. Spectral regions 
desirable for a satellite experiment are not necessarily the ones used here for examples. 
In an actual experiment, two spectral regions could be chosen to measure carbon dioxide 
emission and, . thus, extend the results to higher and lower pressures. A method of deter- 
mining the initial pressure could be developed with two intervals much better than with 
one interval; thus, no extensive work has been done involving only one interval. 

The technique to infer water -vapor structure is applicable to any constituent of 
the atmosphere which has a measurable thermal emission and for which transmittance 
can be calculated. This technique is readily applicable to ozone (O 3 ) which can be deter- 
mined simultaneously with water vapor. 

EXAMPLES OF INFERENCE 

Figure 4 illustrates the inference of temperature from an unperturbed radiance 
profile calculated by using the 1962 U.S. Standard Atmosphere (fig. 2) and temperatures 
inferred from the same radiance profile perturbed with scale, bias, and random radiance 
errors. Pressure ranges of the inferred values are 0.3 to 100 millibars for temperature. 
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Errors used are not unreasonable with present radiometric technology. In the unper- 
turbed inference, the temperature error at the smallest pressure is caused by assuming 
an incorrect lapse rate for the atmosphere above the first level for which temperature 
was inferred. A lapse rate of 4° K/km was used when the lapse rate of the standard 
atmosphere is changing from 2° K/km to 4° K/km. Perturbed examples do not show 
excessive sensitivity to radiance errors. The diameter of the symbol used in figure 4 
is about 2° K. Figure 4(e) illustrates the effect of an error in the assumed initial pres- 
sure P 0 . The error in P 0 of 10 percent produces larger temperature errors than 
the radiance errors shown and indicates that care must be taken in developing a method 
to start the inference. 

Inference of water -vapor mixing ratio is shown in figure 5 for an unperturbed 
radiance profile and for scale, bias, and random radiance errors. Pressure ranges of 
the inferred values are 3 to 200 millibars for water-vapor mixing ratio. Errors intro- 
duced are the same as those used in the temperature examples. Scale error produces 
the largest errors in inferred mixing ratio, and these errors become greater at larger 
pressures. 


CONCLUDING REMARKS 

Techniques to infer stratospheric temperature and water-vapor mixing ratio as a 
function of pressure from limb radiance profiles have been presented. Requirements for 
accurate pointing information have been eliminated, and only the easily obtained vertical 
scan rate is needed. Example calculations of inferred temperatures and water -vapor 
mixing ratio from unperturbed and perturbed radiance profiles are given. Radiance pro- 
files for the 615 cm - -*- to 715 cm“l spectral region due to thermal emission of carbon 
dioxide were used for temperature inference, and radiance profiles for the 205 cm - l to 
295 cm - l spectral region due to thermal emission of water vapor were used for mixing- 
ratio inference. Effects of scale, bias, and random radiance errors are shown. Pres- 
sure ranges of the inferred values are 0.3 to 100 millibars (30 to 10000 N/m2) for tem- 
perature and 3 to 200 millibars for water -vapor mixing ratio. Errors in inferred 
quantities due to radiance errors are not excessive. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., November 25, 1969. 
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Pressure , millibars 



Figure 2,- Comparison of limb radiance profile computed by using altitude with that computed by using pressure 

for 615-cra' 1 to 715-cm~l spectral region. 
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Temperature ,°K 

'(a) No radiance error. 

Figure 4.-' Inferred temperature profiles. 
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(e) Initial pressure error. 
Figure 4.- Concluded. 
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(d) Random radiance error. 
Figure 5.- Concluded. 
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